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Abstract 

We study echoes and what we call 'revival echoes' for a collection of atoms that are described 
by a single quantum wavefunction and are confined in a weakly anharmonic trap. The echoes and 
revival echoes are induced by applying two, successive temporally localized potential perturbations 
to the confining potential, one at time t = 0, and a smaller one at time t = t. Pulse-like responses in 
the expectation value of position {x{t)) are predicted at t ~ nr (n = 2, 3, . . . ) and are particularly 
evident at t ~ 2r. While such echoes are familiar from previous work, a result of our study is the 
finding of 'revival echoes'. Revivals (but not echoes) occur even if the second perturbation is absent. 
In particular, in the absence of the second perturbation, the response to the first perturbation dies 
away, but then reassembles, producing a response at revival times mTx (m = 1,2,...). The 
existence of such revivals is due to the discreteness of the quantum levels in a weakly anharmonic 
potential, and has been well-studied previously. If we now include the second perturbation at 
t = r, we find temporally localized responses, revival echoes, both before and after t ~ mT^, 
e.g., at t ~ mTx — nr (pre-revival echoes) and at t ~ rnTx + nr, (post-revival echoes), where m 
and n are 1,2,... . One notable point is that, depending on the form of the perturbations, the 
'principal' revival echoes at t ~ T^,. it r can be much larger than the echo at t ~ 2r. We develop a 
perturbative model for these phenomena, and compare its predictions to the numerical solutions of 
the time-dependent Schrodinger Equation. The scaling of the size of the various echoes and revival 
echoes as a function of the symmetry of the perturbations applied at t = and t = t, and of the 
size of the external perturbations is investigated. The quantum recurrence and revival echoes are 
also present in higher moments of position, {xP{t)), p > 1. Recurrences are present at t ~ mTx/ j, 
and dominant pre and post-revival echoes occur at fractional shifts of r, i.e. t ~ {mTx i r)/j, 
where the m = 1,2,... and the integer values of j are determined by p. Additionally, we use 
the Gross-Pitaevskii Equation to study the effect of atom-atom interactions on these phenomena. 
We find that echoes and revival echoes become more difficult to discern as the size of the second 
perturbation is increased and/or as the atom-atom interactions become stronger. 

PACS numbers: 03.65.-w, 03.65.Ge 
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I. INTRODUCTION 



In both classical and quantum systems with nonlinearity, a distribution of oscillation 
frequencies leads to dephasing of the response to a temporally localized system perturbation 
and to damping in physical observables of the response. After this response damps away, 
application of another external perturbation can induce subsequent build up of phase co- 



herence, causing later temporally localized responses, ca 
investigated in nuclear magnetic resonance (spin echoes) 



physics js-S], cavity quantum electrodynamics 



m 



led echoes. Examples have been 
li, coupled oscillators J2[j_plasma 



and cold atom systems 



13|. Ad- 



ditionally, in quantum mechanical systems, another phenomena, the so called revivals due 



14|-|l6l|. Revivals are also present 



to the discreetness of the energy eigenstates, can occur 
in systems with a Jaynes-Cummings type interaction, and have been studied both theoret- 



ically 

n 



17 



igj, cavity QED |20|, and circuit QED 



and experimentally in trapped-ion 
|2l| systems. In particular, revivals are macroscopic responses that can result due to re- 
construction of phase coherence in the absence of the second perturbation. Previous work 



has investigated, both experimentally and theoretically, revival [IJ], and echo phenomena 
in trapped atomic systems induced by changes in depth or translations of 10, 11|, the 
confining potential. 

Here we consider echo phenomena induced by two impulsive, successive external pertur- 
bations applied to a collection of atoms in a weakly anharmonic trap. For example, the type 
of external perturbation that we consider in most detail is the application of two translations 
to the trapping potential, one at time t = 0, and a smaller translation at time t = t. We 
study the response of the expectation value of position {x{t)). In the case of a Bose-Einstein 
condensate confined by the pondermotive force of a laser beam, translations of the potential 
can be realized by switching of the laser beam position, or by a change in phase of two 
interfering beams 10| . As in previous work 



ll| , this method results in the creation of 



echoes, e.g., at t ~ 2r, and in the quantum mechanical revival described in |14|-|l6| at, e.g., 
t ^ Tx- In addition, however, we also find quantum 'revival echoes' occurring both before 
and after the revival, e.g., at t ~ T^. ± r and t ^ T^. ± 2r. Furthermore, we find that these 
revival echoes can have much larger amplitudes than, say the echo at t ^ 2r (conditions for 
this to be the case are discussed in Sec. IIIip . In Sec. [TTlwe develop a quantum mechanical, 
perturbative model (in the limit of small anharmonicity and small second displacement) for 



echoes and revival echoes. We then compare these results to the numerically solved, time 
dependent Schrodinger Equation. We study the behavior of echoes and revival echoes at 
different choices of anharmonicity and time delay r, as well as how the size of the various 
echoes scale with the size of the external perturbations. In Sec. IIIII we discuss what effect 
the spatial symmetry of the external stimuli has on the echo phenomena, in particular the 
echoes' relative sizes. In Sec. llVl the quantum recurrence and revival echoes are investigated 
in higher moments of position, {x^{t)), p > 1. Recurrences are present at t ~ mTx/j, and 
dominant pre and post-revival echoes occur at fractional shifts of r, i.e. t ^ {rriTx ± T)/j, 
where m = 1,2,... and the integer values of j are determined by p. Lastly, we model 
atom-atom interactions in our system using the Gross-Pitaevskii Equation, and study the 
effect of interactions on the revival echo phenomena (Sec. IVT) . 

The previous references most closely related to our work are the Refs. js 
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on echoes 

in systems of confined atoms, none of which find our revival echoes. Reference ^ presents 

a theoretical treatment of echoes in confined, cold atoms induced by successive changes in 

the depth of the confining potential. Reference [9| theoretically and numerically studies 

the echo phenomena in a Bose-Einstein condensate where the dephasing is reversed with 

an external optical potential. Reference 10| reports on experimental observations of echoes 

induced by sudden shifts of the potential well realized by changing the relative phasing 

of interfering laser beams, while Ref. Ill] experimentally and numerically compares the 

effectiveness of translations with different temporal profiles in creating echoes. Reference 
□ 

[121 ] experimentally observes an echo in a collection of trapped atoms induced by applying 
a short microwave vr pulse, where the echo is measured via Ramsey spectroscopy. Reference 
[l3 ] numerically implements a Bose-Hubbard model to describe a collection of atoms in an 
optical lattice and their echo response to pulses of radiation in the presence of atom-atom 
interactions. 



II. MODEL AND NUMERICAL RESULTS 

We begin by looking at echo phenomena present in the absence of atom-atom interactions. 
We study the quantum evolution of an initially Gaussian state in a weakly anharmonic, one- 
dimensional trap. 



z^V^ = (ifo + i^i)^, (1) 

Hi = \px\ (3) 

where /3 quantifies the anharmonicity of the trap, with /3 <^ 1. We use harmonic units, 
X = X / {^Jh / mujo) , t = Wot, where x and i are the unnormahzed units, ujq is the frequency 
of the harmonic oscillator, and m is the mass. The procedure will be the following. 

1. Begin with a Gaussian wavepacket centered x = 0. 

2. At t = translate the state by an amount di, x ^ x + di. 

3. At t = T translate the state again in x by an amount d2- 

Later (Sec. IIIip we will discuss another type of perturbation that is not a translation, 
and we will find that the spatial symmetry of the perturbations can have important effects. 
In order to most clearly stimulate echoes, revivals, and revival echoes, we will initially take 
6/2 ^ 1- We first seek the frequencies a;„ of our anharmonic Hamiltonian, H = Hq + Hi, for 
Hi (i.e., (3) small. Expanding Un to second order in /3, we write 

c.„ = a;f + 5c.« + 5a;i^). (4) 

The corrections to the frequencies are conveniently calculated using the creation/ annihila- 
tion formalism of the harmonic oscillator. 



^ = Y^(«^ + «), (5) 
p = ^Y^(«^-«), (6) 

d\n) = y/n\n - 1), (7) 



d'^\n) = Vn + T\n + 1), (8) 



Taking Hi = (/3/16)(a''' + a)^, perturbation theory |22| yields 



= {n\Hi\n) = ^(3 + Qn{l + n)), (9) 



5 



256 



\{n\H, 



(n + 1/2) - (A; + 1/2) 

[2n(l - 2nf{n - 1)] Q{n - 2) + 



-n(n^ — 6n^ + lln — 6) 
4 



e(n-4) (10) 



- 2(2 + 3n + n^){3 + 2nf - -(2 + 3n + n2)(n2 + 7n + I2f \, 



where HQ\n) = ujn^\n) defines the state |n) of the unperturbed harmonic oscillator, and 
0(n) = 1 for n > and Q{n) = for n < 0. The position translations described in the 
procedure can be implemented using the unitary translation operator, 



T{d) = exp[—ipd], 



(11) 



where d is the displacement size, d = di ot d = d2- The matrix elements of this operator in 
the harmonic oscillator basis are 



{■m\T{d)\n) = exp 



11 
' 2 



E 

q=0 



(m — n + g)!(n — q)\q\ 



6(m — n + q), 



(12) 



where 7 = rf/ \/2. There are two special cases of Eq. f|T2|) that will be useful in our analysis. 
The first is the displacement of the harmonic oscillator ground state \n) = |0) (a Gaussian 
wavepacket centered at x = 0), 



(m|T(d)|0) = exp 



7_ 
" 2 



7' 



= exp 



11 
' 2 



C(m,7), 



(13) 



where C(m, 7) = 7™'/vm!. This expression is the well known coefficient for a coherent state 
formed by displacing the harmonic oscillator ground state. 

The second useful expression is the case when 'jy/n <^ 1. In this case, we can approximate 
the matrix element as 



{■m\T{d)\n) = exp 



11 
' 2 



5m,n + 6m,n+llVn+ 1 - 6„i.n-ll VnQ{n - 1) 



Sm,nl'nein - 1) + -^^7V(^ + l)(^ + 2) 



+ ^!^72^/^;r^0(^ _ 2) + (9(7^) 



(14) 
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noting that small displacements predominantly excite nearby states. 

We first consider the response of the system due to an initial displacement at t = 0. 
Beginning with the harmonic oscillator ground state = |0)) at t = 0, we apply a dis- 
placement di. Using the described approximations for the frequencies, and, approximating 
the eigenstates as the unperturbed states \n), for a time < t < r, the quantum state is 



\ip{t < r)) = exp 



y^C(7T,,7i) exp[-iujrit\\n). (15) 



After the displacement di, the probability of being in the n-th state is exp[— 7^]C(n, 71)^, a 
Poisson distribution with mean, 



and variance 



n = Yi, (16) 



< = n = it (17) 



(Physically, the value of n is given as the classical energy over the harmonic oscillator 
frequency, = 1, ri = (p^/2 + /2)/ujq = d\/2 = 7^.) 

The response of the system, as measured by the expectation value, {ilj\x\ip), is given by 



C{n)C{n + l)Vn+ lcos[(a;„+i - a;„)t], (18) 



where for notational clarity we have dropped 71 in the argument of C{n) = C(n,7i), and 
A = exp[-^f/2]. 

For comparison, we obtained numerical solutions of the Schrodinger Equation using a 



split-step operator method (23|]. Figure [It^a) shows a plot of {x{t)) versus t obtained by 
numerical solution of the Schrodinger Equation in the absence of the second displacement. 
We see from Fig. [T](a) that the displacement at t = leads to a fast oscillatory response 
whose envelope decays to small values by t ~ 1500. At longer time, the response reassembles 
with its envelope reaching a peak value at the 'revival time' = 8625. This revival of the 
response has been called a quantum revival or recurrence. (Note that in Fig. [T] the fast 
oscillations are not visible because the thickness of the line that is plotted exceeds the 
period of the fast oscillations.) Figure (H^b) plots the value of {x(t)) calculated from Eq. 
(fT8|) with Eqs. (j4]), (Q and (fTOj) used for the frequencies. We note very good agreement 
between the model and the numerical solution. 
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14-16 



24 



. (The 



25j. This 



As we now verify, this term includes the well known quantum recurrence 
expression for Eq. f|T8|) does not capture the fractional revival described in 
revival is much smaller in amplitude than the large revival at t ~ T^., and is not visible on 
the scale used in Fig. [T] We have verified that a model including 0{P) corrections to the 
energy eigenstates does reproduce these small fractional revivals). In order to estimate the 
quantum recurrence time of {x(t)), in a manner similar to Refs. 15|, ll6|, we expand our 
expression for the frequency difference around a mean value of n, denoted by n. Using Eqs. 
dl]) and (fTOj) . the energy difference between two adjacent states can be approximated as 

Un+i - ujn ^ ^uJn + a{n - n) + (19) 

where 

Ac. = 1 + l^in + 1) + + 1^ + 1^') (20) 

a = ^(8/3-17/32-17n/32). (21) 

Aufi is the difference in the two energies, evaluated at n (and is independent of n). We can 
then express cos[(co'„+i — C0n)t] as: 

cos[(Aa;fj + a{n — n))t] = Re{ exp[iAci;fit] exp[ia(?T. — n))t]Y (22) 

In the vicinity of n, the oscillatory behavior is a product of a fast oscillation at a frequency 
Aufi, which is independent of n, and a slow oscillation, which is periodic for all n at a period 
27i/a. Thus the envelopes are periodic at times t = mTx {m = 1, 2, ...) where 

27r 



a 



(23) 



which agrees with the value = 8625 estimated by inspection of Fig. [TJ^a). 

It should also be noted that the time ~ T^ev/'^i where T^e,, is the revival time it 
takes the state to come back to approximately its initial condition. The recurrence in (x(t)) 
corresponds to a mirror revival [16|, e.g.: the quantum state has reassembled as a mirror 
image of its initial condition on the opposite side of the anharmonic well. 

Now, including the effect of the second displacement, the state is allowed to evolve to time 
t = T, at which time it experiences a small second displacement by an amount ^2. In order 
to utilize the expression in Eq. ( IT^ we define 72 = d^j^^ and require 72 ^ 1- Taking as 

8 
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t 



FIG. 1. (a): {x{t)) found through numerical solution of the Schrodinger Equation for a system 
with only one initial displacement at t = 0, /3 = 0.001, di=5. (b): {x{t)) found through the 
perturbative model, for the same parameters as (a). Note good agreement between the model and 
the numerical solution. 



an upper bound on n to be n + 3cr„,, and noting that n = 7^ we require 72 <^ (71 + 871) 
Using the expression for {'m\T{d)\n) in Eq. ([T^ we find 

9 



\ip{t > r)) = exp 



7i 
' 2 



'^^{m\T{d2)\n)C{n) exp[— zcunr] exp[—iuJmit — 'r)]\m) 

m,n 

A\ C{m) exp(— icjm^)!^) 

m=0 

+ 72 ^ C'l"^ — l)A/mexp[— zcjm,_ir] exp[— zw^li — 

m=l 

oo 

— 72 C(m + l)v^m~TTexp[— zWrn+iT] exp[— — 

m=0 

oo 

— 72 C(m)mexp[— zwrn^ll^) 



m=l 



12 



C{m — 2)^/m{m — 1) exp[—i(jJm~2T] Gxp[—iUm(t — r)]\m) 



m=2 

oo 



2 °° 1 

+ y XI ^^"^ + 2) V(m + l)(m + 2) exp[-iw^+2r] expHu;„(t - r)]|m) + 0(72) L 

m=0 

(24) 



where A = exp [—(7^ + 7|)/2]. Using the expression for the quantum state for t > r, we 
can calculate the expectation value of the position {x{t)) = {tp\x\ip) and order the terms in 
according to their dependence on the second displacement 72, 



(1) 



(25) 



The 0th order term {x{t))^^^ is the result in the absence of the second displacement and is 
given by Eq. ( IT8l) . Calculating the next highest term in 72 in our sum, we find an expression 
for our first revival echo, which can be separated into two parts, 



(x(t))W = {xY^' + {xy_^i 



(1) , 



(26) 



where 



1 jx)? 
72 AyV2 



2 X(n + 1) {C{nf - C{n + 1)^} cos[(a;„+i - a;„)(^ " r)], (27) 



n=0 
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and 



1 / °° 



cos 



— COS 



[Un+l - UJn)t + {Un - COn-l)T 

(28) 



The superscript refers to the power of 72 (ignoring the common dependence in the nor- 
mahzation constant A), and the subscript indicates the time of the echo, in units of r, 
relative to t = T^. For example, the first sum has a cosine dependence similar to the 

quantum recurrence in Eq. ( ITSl) . except t has been replaced with t — t. Thus, we expect 
a large amplitude in this second term at t = r and t mT^ + t (m > 1), which we call 
a post-revival echo. The second sum (x)*^] has a time dependence which again depends on 
the difference of adjacent frequencies. Note however, that if we take the spacing between 
adjacent energy levels to be approximately independent of n, Un+i — uJn ^ Un — Wn-i, then 
cos[(ci;„+i — uJn)t + {oJn — Ci;„_i)r] ~ cos[(ci;„+i — uJn){t + r)]. Thus, we expect the second 
sum to contribute to a pre-revival echo at t ~ TnTa^ — r. Note that our expression for 
only applies for t > t. Thus the m and r values for pre-revival echoes are restricted by the 
requirement that mTj, > 2r (e.g., for r < we have that m > 1 is permissible). 

Continuing on to {x{t))^'^\ we can again separate this term by the location of the dominant 
echoes, 

{x{t)y'^ = (^)?^ + (^)f + (29) 

where expressions for {x)''q \ {x)^\ and (x)^2 given in Appendix Rl 

With a similar argument for the time dependence of the first order terms, approximating 

the spacing between levels to be almost constant, each sum in (x)q produces a second order 

(2) 

response at t ^ T^. Likewise, the second term (x)2 produces an echo at t ^ 2r and and a 
post-revival echo at t ~ rriTx + 2r, while pre-revival echoes from the last term (x)_2 occur 
at t ^ mTx — 2r (m > 1). 

Figures [2] and [3] compare numerical solutions of the Schrodinger Equation (Figs. [2]J^a,b) 
and Fig. [3](a)) with predictions (Figs. [2](c,d) and Fig. Mjo)) of the perturbation theory 
model [i.e., Eqs. ([H]), ([26]) -([28]), (jATll-dAll) for {x{t)), with Eqs. (HD, and ^ for u^] 

11 



for different choices in /3, r, and d2 parameters. In Fig. [21 is calculated for /3 = 0.001, 

T = 1499, d2 = 0.05. There is good agreement between the numerics and the model, with 
the model reconstructing the quantum recurrence at t ~ Tj., and the echo and revival echoes 
at t ~ 2r, t Tx — T and t ^ — 2r. 

Similarly, Fig. [3] plots {x{t)) for a larger /3 = 0.002 and d2 = 0.1, and a smaller value 
of r = 1200, showing more than one quantum revival echo. There is again good agreement 
between the model and the numerics in capturing the behavior of the quantum recurrence 
at t ^ mTx {m = 1,2) and the echo and revival echoes at t ^ 2r, t ^ + r, t ~ mTx — r 
and t ~ mTri: — 2r (m = 1, 2). 

Additionally, the model indicates that one should expect the amplitude of the pre-revival 
echo a t ~ — r to scale linearly with the second displacement d2- As seen in Eq. ( 1271) the 
term (x)^^"* is proportional to 72 (and therefore ^2) if we ignore the weak dependence due 
to the normalization constant A. Similarly, the echo at t ~ 2r and the pre-revival echo at 
t ^ Tx — 2t (represented by (x)2^^ and {x)^^l) are expected to increase quadratically with ^2, 
as shown in Eqs. flASP and (IA4p of Appendix |X1 We study the scaling of these echoes using 
the full numerical Schrodinger Equation with the choice oi (3 = 0.001, r = 1899, di = 5, and 
various values of d2. This value of r is chosen to avoid the coincidence of the t ^ 2t echo 



and the t ^ — 2t pre-revival echo with the small fractional revival described in 2S, |25 |. 
The amplitudes of the echoes and pre-revival echoes as a function of ^2 are plotted in Fig. HJ 
and compared to the amplitude scaling predictions of the perturbative model. We note that 
indeed the pre-revival echo at t ~ T^. — r appears to behave linearly with ^2, while the echo 
at t ~ 2r and the pre-revival echo t ^ — 2t appear to grow quadratically, as expected. 
In general, the model suggests that echoes at t ~ nr and revival echoes at t ^ mTx ± nr 
(m > 1, n > 1) should have lowest order scalings of ((^2)". (It should be noted, however, 
that the scaling of very small echoes may be affected by revivals and echoes occurring due 
to corrections in the energy eigenstates.) 

An approximate dependence of the pre-revival echo at t ~ T^; — r on the first perturbation 
di can also be obtained, and can be understood by studying Eq. ( 1271) . At the time of maximal 
{x(t)) (in the vicinity of t ~ T^, — r), one expects that for n near n, the cosine terms will be 
approximately in phase and only depend weakly on n. Thus, taking the cosine terms to be 

12 




FIG. 2. (Color Online): (a): {x{t)) solved numerically for the full Schrodinger Equation for 
/3 = 0.001, di = 5, ^2 = 0.05, r = 1499. (b): same as (a), but zoomed in to display the smaller 
echoes, (c): The model {x{t)) plotted for the same parameters, (d): same as (c), but zoomed in 
to display the smaller echoes. We note that the model reconstructs the t ~Tx — t pre-revival echo, 
as well as the echo at t ~ 2t and the pre-revival echo at t ~ — 2r. 

roughly constant, and independent of n, 

oo 

{xt\ ~ 72 exp[-7?] eM-ll] E + " + 1) (30) 

n=l 

and using definitions of C(n), A, and exp[7^] = X]^o(7i)"/''^'' ^"^^ 

~ 7?72 exp[-7^] ~ 7^72 ~ d\d2. (31) 
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FIG. 3. (Color Online): (a) {x{t)) solved numerically for the full Schrodinger Equation for 
P = 0.002, di = 5, ^2 = 0.1, T = 1200. (b) The model {x{t)) plotted for the same parameters. 
The model reconstructs the quantum recurrence at t ~ rriTx, as well as the echo and revival echo 
responses at t ^ 2t, t^Tx + T,t^Tx — t, and t ^ Tx — 2t, (m = 1, 2). 



Note that we ignore any dependence on di contained in the cosine terms (arising from 
dependence on n of the expansion of Un+i — u^n, which in turn is related to di). Figure [5] 
plots the amplitude of the pre-revival echo as a function of dld2. Note that echo amplitude 
appears to scale linearly with dld2. For comparison, in the inset we plot the amplitude as a 
function of ^2, where lines connect simulations with equal values of di. 

A similar procedure can be implemented for the response at t ~ r, and the echoes at 
t ^ Tx + T, t ^ 2r and t ^ T^. ± 2r. The response at t ~ r and the revival echo at t ~ T^, + r 
are given by Eq. fl27|) . Summing over the coefficient gives 

{x)?-d2 (32) 

with approximately no di dependence, which we numerically verify but do not show. For 
the echo at t ^ 2t and the revival echoes at t T^. ± 2r, the oscillatory terms in Eqs. f lASP 



and (IA4p are assumed to be in phase, and the summations are done over the coefficients, 
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-1 



FIG. 4. (Color Online): The amplitude of three echoes present at different times as a function of d2 
for (ii = 5, /3 = 0.001, r = 1899, {t ~ — r: circles, t k,Tx — It: squares, t ~ 2t: triangles). The 
solid lines are echo response amplitudes extracted from the perturbation theory model. The dotted 
lines are guides to the eye indicating linear and quadratic behavior. We note that the pre-revival 
echo at t ~ Ta; — r scales linearly with d2 , while the echo responses at t ~ 2t and t ^ — 2t scale 
quadratically. 



Thus we expect the echoes and revival echoes to approximately scale as (i^c?!- Figure El 
which plots the echo amplitudes as functions of clfd"^, for /3 = 0.001, 3.5 < di < 6 and 
0.001 < d2 < 0.07 and r = 1699, confirms this expectation. (The post-revival echo at 
t ^ Tj, + 2r is the recurrence of the echo at t ^ 2r. It has an amplitude very similar to the 
t ^ 2t echo and is not shown in Fig. [6] for clarity.) 

Increasing d2 to be ~ 1 (where we no longer expect our perturbative model to be valid), 
Fig. [7] plots the values of {x{t)) for /3 = 0.001, r = 1499, di = 5, and different values of 
d2 ranging from ^2 = 0.1 to ^2 = 0.5. As d2 is increased, we see that the amplitude of the 



leading to 



^ ~ dfdj + 0{didl) 



(33) 
(34) 
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FIG. 5. (Color Online): The amplitude of the pre-revival echo at t w T^,. — r as a function of (ifd2 
for /3 = 0.001, r = 2499, and various values of di {di = 3.5: circles, di = 4.0: squares, di = 4.5: 
asterisks, di = 5.0: triangles, di = 5.5: stars, di = 6.0: x's). d2 is allowed to vary between 
0.001 < d2 < 0.07. The dotted lines are echo amplitudes extracted from the model equations. We 
note that the pre-revival echo at i ~ T^; — r scales linearly with did2. (Inset): The amplitude of 
the revival echo as a function of d2, where lines link simulations done at equal values of di. 

quantum recurrence at t ~ Tj, diminishes. Further, the amplitude of the pre-revival echo at 
t ^ Tx — T is also suppressed, while the amplitude of (x(t)) increases for shorter times. 

Finally, we note that from numerical experiments solving the Schrodinger Equation, at 
larger di we clearly see an echo at t ~ 3r and revival echo at t ~ — 3r. An example 
illustrating this is shown in Fig. [HI 
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FIG. 6. (Color Online): The amplitude of the echo at t ~ 2r (red, solid) and the pre-revival echo 
at t ~ Tx — 2t (blue, unfilled) as a function of (i'^dg for, /3 = 0.001, r = 1699, and various values of 
di and d2- 3.5 < di < 6 and 0.001 < fi2 < 0.07. The dotted lines are echo amplitudes extracted 
from the model equations. We note that the t ~ 2r echo and the pre-revival echo at t ^ Tx — 2t 
scale linearly with dfdl- 
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FIG. 7. (Color Online): {x{t)) for /3 = 0.001, r = 1499, di = 5, and different values of ^2, (a) 
d2 = 0.1, (b) d2 = 0.2, (c) d2 = 0.3, (d) da = 0.4, and (e) = 0.5. Note that both the 
quantum recurrence at t ~ T^, and the pre-revival echo at i ~ — r are suppressed as da increases. 
Further, the amplitude of {x{t)) is larger for shorter times as da increases. 
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FIG. 8. (Color Online): (a) and (b): {x{t)) sovled numerically for the full Schrodinger Equation 
for /3 = 0.001, di = 8, d2 = 0.05 and r = 1299. Note that an echo at t ~ 3r and pre-revival echo 
t K Tx — 3t are clearly evident, (b) is the same as (a), plotted on a different scale to better display 
the echoes. 
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III. DEPENDENCE ON THE SYMMETRY OF THE EXTERNAL PERTURBA- 
TIONS 



In Sec. ini we considered echoes and revival echoes induced by two successive displace- 
ments of a the trapping potential by amounts di and d2- For ^2 small, we found that the 
revival echoes at t ~ mTx ± r can be much larger (i.e., 0{d2) as opposed to (9(^2)) than 
the first ordinary echo (i.e., the echo occurring at t ~ 2r). A possibly significant point is 
that for small displacement d and symmetric traps, V{x) = V{—x), the perturbation may 
be viewed as being antisymmetric: V{x + d) — V{x) = d{dV/dx) + O^d"^). Thus, a natural 
question is whether our result for the relative sizes of the t ^ 2t echo and the revival echoes 
at t ~ TnT^ ± r depend on the symmetry of the stimuli. In what follows we consider this 
question and show that the answer is affirmative. 

We have considered a different type of perturbation at t = r which is symmetric in x, 
an impulse squeeze, i.e., we add a term a2X^6{t — r) to the potential. Taking /3 = 0.001 
and T = 1299, we numerically calculate the response to an initial shift {di = 5) and a small 
impulse squeeze (0^2 = 0.005). We display the results in Fig. |9l In this case, echoes and 
revival echoes are evident in {x{t)), at t ~ 2r and t ^ T^. — 2r, with no responses in {x(t)) at 
t = T and t ^ Tx~T. Model expressions for the echoes can be found in the same manner as 
in Sec. [Tll and predict that the t ^ 2t echo and t ^ — 2r revival echo both scale linearly 
with 0^2 (see Appendix [B|) . Figure [TO] plots the amplitude of the echo as a function of ^2, 
demonstrating that the echo amplitudes depend linearly on 02- This behavior is due to the 
impulse squeeze only exciting states of like parity; e.g., the impulse squeeze acting on an 
even eigenstate only excites the other even eigenstates. Note that unlike our result in Sec. 
IIT| the first {x{t)) echo (i.e., at t ~ 2r) and the first revival echo (i.e., that at t ^ — 2r) 
both scale linearly with the strength 02 of the second perturbation. (In fact it can be shown 
that all revival echoes at t ~ friT^ ± 2r scale linearly with 02.) 
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FIG. 9. (Color Online): {x{t)) solved numerically for the full Schrodinger Equation with a dis- 
placement at t=0 (di = 5), and an impulse squeeze at t = r = 1299 (02 = 0.005). Notice there is 
no response at tr or a pre-revival echo at t ^ — t. An echo at t ~ 2r and a pre-revival echo at 
t ~ Ta; — 2t are evident and are comparable in size. 
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FIG. 10. (Color Online): The amplitude of the echo at t ~ 2r (triangles) and the pre-revival echo 
at t w Tx — It (squares) as a function of 02 for, (ii = 5, /3 = 0.001, r = 1299, and various values 
of 02. The dotted line is a guide to the eye indicating linear behavior. We note that both echoes 
scale linearly with 02. The solid lines are echo amplitudes extracted from a perturbation theory 
model (see Appendix [B]) . 
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IV. QUANTUM RECURRENCES AND REVIVAL ECHOES IN (x^ (t)) 



Returning to the response to two successive displacements (as in Sec. [TTl), we now examine 
the behavior of the revival and the revival echoes in the response of {x^{t)), for p > 1. 
For example, if one were to measure the width of the cloud of cold atoms, the necessary 
observable would be w'^{t) = (x^(t)) — Performing an analysis analogous to that 

of Sec. ini (described in detail in Appendix [C]) , we find that quantum recurrences (in the 
absence of the second perturbation) are to be expected at fractional times oi t = T^. In 
particular, for the observable {x^{t)) , recurrences are expected at times t ~ mTx/j with 
m > 1. The values of j are contingent on whether p is even or odd; for p even: j = 2,4, . . . ,p, 
and for p odd: j = 1,3, ... ,p. The amplitude of the recurrences is smaller for larger values 
of j. 

Additionally, revival echoes are also present in (x^(t)). The dominant echoes (those that 
scale as ^2), are also described by our model in Appendix O Pre-revival echoes are expected 
at t ~ {friTx — t)/] and post-revivals at t ~ {fnTx + T)/j, where again m > 1 and the 
possible values of j are different for even and odd p; for p even: j = 2,4, ... ,p, and for p 
odd: j = 1,3, . . . ,p. For example, for p = 3, recurrences are expected at t ~ mTx/3, as well 
as larger revivals at t ~ TnT^, m > 1. Additionally there are dominant pre and post-revival 
echoes at t ~ {fnTx ± r)/3 and t ^ mT^ ± r, for m > 1 . The model suggests that echoes 
at t ~ nr/j, (t > r), and revival echoes at t ~ {fnTx ± nr)/ j will scale as 72 . 

Figures fTTT a) and fT2r a) display the values of {x^{t)) calculated by solving the Schrodinger 
Equation, for p = 2, and p = 3 after an initial shift perturbation, di = 5, at t = 0, and a 
second shift perturbation, d2 = 0.05, at t = r = 1499, with /3 = 0.001. Displayed in solid 
lines are the expected times of the quantum recurrences, and in dashed lines, the expected 
times of the pre and post-revival echoes. For p = 2, echoes at t ^ 3r/2 and t ^ 4r/2 
are expected to be very small, scaling as (^2)'^ and (^2)^, and are not readily evident upon 
examination of Figure [TTT a). Regarding the scaling, as for the p = 1 case, we note that 
the scaling of very small echoes may be affected by revivals and echoes occurring due to 
corrections in the energy eigenstates. 

In Fig [T2I the pre-revival at t ^ T^- — r/3 is not visible because of the large quantum 
recurrence at t ~ T^. A conventional echo is expected at t ~ 2r, scaling as (^2)^ but is 
obscured by the revival at t ~ Tx/3. For larger values of d2, and shorter values of r, we have 
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numerically observed the t ^ 2t echo. 

Using the model expressions described in Appendix [Cl we compare the numerically cal- 
culated values of (a;^(t)) and {x^{t)) to those obtained from our model, as shown in Figs. 
[TTjfb) and [T^ b). Note that there is good agreement between the numerics and the model. 
In both the {x^{t)) and {x^{t)) cases, the response is dominated by the quantum recurrences, 
and the pre and post-revival echoes. 



T /2-T/2 T /2 T /2+T/2 T -t/2 T 

XXX XX 




FIG. 11. (Color Online): (a) and (b): sovled numerically for the full Schrodinger 

Equation for /? = 0.001, di = 5, d2 = 0.05 and r = 1499. Note that recurrences are present at 
t ^ Tx and t ~ Tx/2 (solid lines). Also, pre-revival echoes at t ~ {mTx — t)/2, m = 1,2 and a post 
revival echo t ~ {Tx + t)/2 are evident (dotted lines), (b) The value of {x'^{t)) calculated using the 
model in Appendix [Cl 
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FIG. 12. (Color Online): (a) and (b): {x^{t)) sovled numerically for the full Schrodinger 
Equation for f3 = 0.001, di = 5, ^2 = 0.05 and r = 1499. Note that recurrences are present at t ~ T^; 
and t ~ mTx/3, m = 1,2 (solid lines). Also, pre-revival echoes can be seen at t ~ (mTx — t)/3, 
m = 1,2 as well as t ~ Tx — t. Post revival echoes are evident at t ~ {mTx +t)/3 , m = 1,2 
(dotted lines), (b) The value of {x^{t)) calculated using the model in Appendix [Cl 



25 



V. ATOM- ATOM INTERACTIONS 



Finally, we include a cubic nonlinear term in the Hamiltonian and numerically solve for 
the wavefunction. This can be considered as the mean field approximation for the condensate 
wavefunction of a Bose-Einstein condensate with interactions, and is known as the Gross- 



Here u is a nonlinear interaction strength, which quantifies repulsive {u > 0) or attractive 
{u < 0) interactions in the condensate. Again considering echo responses to two successive 
displacement (as in Sec. |Tl]) and taking di = 5,d2 = 0.05, r = 2499, and (3 = 0.001 for 
various values of u, we numerically investigate the behavior of the echoes and revival echoes. 
The results are plotted in Fig. [T31 As the interaction strength is increased, we note that the 
damping immediately after t = is increased, as compared to the u = case. Additionally, 
for larger values of u, the value of does not fully damp away before the application 

of the second displacement. This observation has been made previously for similar systems 
29|. The inclusion of interactions also suppresses the quantum recurrence at t ~ Tj., as 
well as the amplitude of the revival echo at t ~ — r. The time at which the quantum 
recurrence occurs is also shifted from the predicted value of t ~ T^. in the noninteracting 
case. 



Pitaevskii Equation (GPE) |26|-l28|. 




(35) 
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FIG. 13. (Color Online): The value of {x{t)) as a function of time with di = 5, d2 = 0.05, /3 = 0.001, 
r = 2499, for various values of u. (a) u = 0.05, (b) u = 0.10, (c) u = 0.15, (d) u = 0.20. We note 
that while increasing interactions causes damping to occur faster, it also leads to nonzero values 
of {x{t)) before the second kick at t = r. Further, as u is increased, both the quantum recurrence 
t ^ Tx and the pre-revival echo at t ^ — t are suppressed. 
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VI. SUMMARY AND CONCLUSIONS 



In this paper, we study echoes and revival echoes for a collection of cold atoms in a 
weakly anharmonic potential subjected to two external stimuli, one at t = and the other 
at t = r. In the case where the two external stimuli are sudden displacements of sizes di at 
t = and at t = r, we observe responses in the expected value of position and find 

that responses occur not only at t ^ nr, {n>\) and at well-known quantum recurrences at 
t rriTrc, but also at t ~ rriTx — nT>T (pre- revival echoes) and t ^ mT^. -|- nr (post-revival 
echoes). Again, in the case where the external stimuli are sudden displacements, we have 
formulated a perturbation theory model and note good agreement between our model and 
the numerical results. We numerically verify that for small ^2, the revival echo at t ~ Tc ±r 
scale linearly with d2 while echoes at t ~ 2r and revival echoes at t ~ T^; ±2r scale as ((^2)^. 
As expected from our perturbative model, and seen numerically, this scaling suggests that, 
for a sufficiently small d2, the revival echoes at t ^ ± r can be substantially larger than 
the echo at t ~ 2r. In addition, our perturbation theory model shows that for small di and 
d2, the revival echo at t ^ — r scales as d1d2, (which we numerically verify), while the 
echo at t ~ 2r and the revival echoes at t ~ Tj. ± 2r scale as d^d2 (which we numerically 
verify) . 

We numerically demonstrate the suppression of the quantum revival t ^ and the 
revival echo at t ^ Tj. — r, when the size of the second displacement is increased, ^2 ~ 1. 
We also investigate how the echoes and revival echoes depend on the form and symmetry of 
the external stimuli. One result is that, in the case where the first stimuli is a displacement, 
the domination of the size of the revival echoes at t ^ TnT^ ± r over the echoes at t ^ 2r 
does not occur when the small second external stimulus has different symmetry than that 
of a small displacement. 

The presence of quantum recurrences and revival echoes in responses observed in higher 
moments of position was also studied. Quantum recurrences of (x^(t)) appear at fractions 
of the p = 1 revival time, t ^ mTx/j, where m = 1,2,..., and the value of j is determined 
by p. Revival echoes are also present in x^, with dominant echoes (scaling as ^2) at {x^{t)) 
at times t ~ {mTx ±T)/j. 

Finally, the the inclusion of interactions, modeled with the Gross-Pitaevskii equation, 
demonstrates that the quantum recurrence and the revival echoes are suppressed as the 
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nonlinear interaction strength u, is increased. 

Imaging is the natural way to resolve revivals of x. Revivals and echoes in momentum, 
p, may be found by a similar calculation. In this case, these revivals and echoes may be 
observed if the trap is switched off at a revival or echo time. The cloud of atoms will move 
together since all atoms have nearly the same momentum. If instead, the trap is turned off 
at a time not coinciding with a revival or echo, the cloud of atoms will simply spread. 
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Appendix A: Expressions for {x)^q \ {x)'^\ and {x)^^l 

Continuing the procedure outlined in Sec. [ITl we can find the terms in the expectation 
value of {x{t)) which are quadratic in 72. These can again be separated into the location of 
their dominant echoes. 



(Al) 



The expressions for for (x)q , {x) 



and {xfl 



are found to be: 



00 



2^C(n + l)(:7(n + 2)(n + l)Vn + 2cos {u, 



'n- 



n=0 



(A2) 



00 



+ 2'^C{n-l)C{n){n + l ) \fn cos [uj, 



'.„+i - a;„)(t - r) + (w„ - a;„-_i)r 



n=l 




n=0 



2 l)cos (, 



n=2 
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J:£Z2 =_2Y C{n)C{n + l){n 



+ 1)^/^ cos - Un){t - r) + - a;„+i)r 



(A3) 



^C(ra)C(n - l)(n + l)v/ricos (wn+i - + - w„+i)r 

n=l 

oo 

C{n)C{n + l)nA/?^ + 1 cos (a;.„ - Wn-i)^ + (w„_i - cj„+i)r 



n=l 



and 

7i 



/^\(2) °° ^ 

= -2 E - 2)C'(^ + l)Vn(n2-l)cos 



(A4) 

oo 

+ ^ C(n)C(n - 3)A/(n)(n - l)(n - 2) cos (w„ - Un^i)t + - Wn-s)^" 



n=3 

oo 



^C(n)C(ra + 3)A/(r2 + l){n + 2)(r2 + 3) cos {un+i - uJn)t + (a;„+3 - cj.„+i)r 



n=0 



In order to get the approximate scaling for the echo at t ~ 2r and the revival echo at 
t ^ Tx — T, we take the above expressions for the echoes, and take the cosine terms to be in 
phase and independent of n. Doing the summations for the coefficients, we find for (a;) 2 • 



C{n)C{n + l)(n + if' = e^llM + 7i) 

n=0 

00 

Y,C{n)C{n - l){n + l)^ = exp[7?](7? + 271; 

n=l 

00 

^C(n)C(n+ l)nv/^^ = exp[72](73). 



(A5) 
(A6) 
(A7) 



n=l 



and for {x)^^l: 



J2 C{n - 2)C{n + l)]^^^^ - 1) = exp[7,2](7D 

n=2 
00 

J2C{n)C{n - 3)v/n(n-l)(n-2) = exp[72](73) 

n=3 

00 

J2 C{n)C{n + 3)nv/(n + l)(n + 2)(n + 3) = exp[72](73). 



(A8) 
(A9) 
(AlO) 



n=0 



The factor exp[7^] is canceled by the 71 dependence of A^, leading to the di dependence 
described in Sec. [Ill 
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Appendix B: Model Expressions for an initial displacement, followed by an impulse 
squeeze. 



In a manner analogous to the calculation presented in Sec. [TTl The response of {x(t)) 
to an initial displacement at t = (rf = rfi) and an impulse squeeze at time t {a = 02 , 
a2{n + 3(T„) ^ 1) can also be estimated. Approximating the matrix elements of the squeeze 
operator in the unperturbed anharmonic oscillator basis, we find 

(m| exp[-ia2X^\n) ^ 6m,n - «Y|(2ra + l)6m,n + ^/n{n - l)6m,n-2 (Bl) 

From this, the responses proportional to a2 can be calculated, 

= {x)S^ + (B2) 

where 

: ;7 — = V C{n)C{n + l)nVn + lsin[(a;„ - a;„_i)t + - UJn+i)^] (B3) 

exp[-7f]a2/v2 ^ 

00 

+ ^ C(n)C(n - + l)v/n sin[(a;„ - cjn+i)^ + (wn+i - Wn-i)r] 



and 



n=l 



(B4) 



r 21 -7^ = E C{n)C{n + 3) ^(n + l)(n + 2)(n + 3) 
exp[-7f]«2/V2 ^ 

X sin[(w„, - w„+i)t + - ujn+s)'^] 

00 

+ EC(n)C(n - 3)^n{n-l){n-2) 

n=3 

X sin[(w„, - ujn^i)t + - a;„_3)r]. 

Appendix C: Model Expressions for {xP{t)) after two shift perturbations. 

In the case of a shift perturbation, (ii, at t = 0, and a smaller shift perturbation, d2, at 
time t = r, we derive model equations for both the quantum recurrence, as well as revival 
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echoes in the quantity 

(x^it)) = ii^m) (CI) 

for p > 1. Using the definition of x in Eq. [5l one can write the matrix elements of x^, in 
the unperturbed harmonic oscillator basis, as 



{m\xP\n) = {^) {m\{a + a^y\n) (C2) 



V2 

1 Y ^ 

-/k) ^B{n,p-2j,p)6m,n+p-2j (C3) 

where B{n,p — 2j,p) are coefficients that can be found via application of the annihilation 
and creation operators. For example, for p = 3, the coefficients are 



B{n, -3, 3) = ^/n{n-l){n-2) (C4) 
B{n, -1,3) =3ny/E (C5) 



B{n,l,3) = 3{n + l)Vn + l (C6) 
B{n, 3, 3) = ^{n + l){n + 2){n + 3). (C7) 

It can be shown that these coefficients obey the relation B{n, —d,p) = B{n — d,d,p). 
Using Eq. (1C3I) . and the expression for the quantum state after two displacements, Eq. 
(^^, one can again find the time dependence of the expectation value of and order 

the terms according to their dependence on 72, 

{xP{t)) = + + . . . (C8) 

The term (x^(t))''°'*is again the quantum recurrence in the absence of the second displace- 
ment, and if p is even, is given by. 



■J 



Un)t]B{n,j,p) (C9) 



J even 

and if p is odd. 



^ ' * ' 1=1 n=0 



j odd 
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Employing the same argument for the behavior of the quantum recurrence in Sec. HIl one 
finds 

^ jAf, + ja(n -n + q). 

where g is a constant independent of n. Compared to the J = 1 case of Sec. [Til the slow 
envelope maxima are now periodic for all n at 27r/(ja). Thus, each cos[(a;„+j — will 
give rise to a recurrence at time 

TflT 

^, (j>l,m>l) (C12) 

If p is even, then the j = case is a constant which increases the value of Note that 

the expectation value of (x^) contains recurrences not only at t ^ rriT^/p, but also at the 
values of j < p as seen in Eqs. (]C9P and (IClOp . 

The expression for the revival echoes that are approximately linear to d2 can also be 
found, for p even: 

and p odd: 

j odd 



j=2 
j even 



^''"^'f'^ = T \{xW-i + i^m?} (C14) 



where, 

')o'^ = $^^(^ + l)C{n)V^^cos[{un+i - u;n)r][B{n + 1, 0,p) - B{n, 0,p)], (C15) 



00 



n=0 



00 



= X]'^*^'^ ~" - jOv^^I'^, -j,p)cos[{ujn-j - uJn)t + (w„ - w„_i)r] (C16) 

00 

^C(n + l)C(r2 + + lB{n,j,p) cos[(a;„+j - + (w„ - a;„+i)r]. 



n=0 

and 



= X]C'('^ - l)C(n + j)v^S(n, cos[K+, - w„)t + K - ^n-i)r] (C17) 

n=l 
00 

^C(n+ l)C(n- + lB{n,-j,p) cos[(a;„_j - ujrn)t + (un - Un+i)r]. 



n=j 
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The term {x^)^q^ is a time independent offset. The quantity represents a post- 

revival echo at time t ^ {mTx + T)/j. Similarly the term gives rise to a pre-revival 

echo at t ~ {mT^ ~ '^)/3 i 3 ^ 1 ^"^^ m > 1. As can be seen from Eqs. flC13l) and 
(lC14p . revival echoes will be present and at various fractional shifts of r from the different 
recurrences. For example, for p = 3, recurrences are expected at t ~ ^^y^, as well as larger 
recurrences at t ^ TnTx-, m > 1. Additionally there are dominant pre and post-revival echoes 
at t ~ {fnTx ± r)/3 and t ^ mTj. ± r, for m > 1 (see Figure [T2|) . 
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